In-Plane Elliptic Flow of Resonance Particles 
in Relativistic Heavy-Ion Collisions 
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We analyze the second Fourier coefficient V2 of the pion azimuthal distribution in non-central 
heavy-ion collisions in a relativistic hydrodynamic model. The exact treatment of the decay kine- 
matics of resonances leads to almost vanishing azimuthal anisotropy of pious near the midrapidity, 
while the matter elliptic flow is in-plane at freeze-out. In addition, we reproduce the rapidity de- 
pendence of for pions measured in non-central Pb -I- Pb collisions at 158^4 GeV. This suggests 
that resonance particles as well as stable particles constitute the in-plane flow and are important 
ingredients for the understanding of the observed pion flow. 

25.75.Ld, 24.10.Nz 



^ The main goals of relativistic nuclear collisions are to determine the nuclear equation of state (EOS) under extreme 



conditions and to understand a new phase of deconfined nuclear matter, the quark-giuon plasma (QGP) Since 
it is the pressure gradient perpendicular to the collision axis that causes various transverse collective flows, such as 
[ — . ■ radial flow, directed flow, and elliptic flow, in relativistic nuclear collisions, these flows observed in the final state are 
, expected to carry the information about the EOS Q. If the QGP phase is created in nuclear collisions, the quark 
matter expands, cools down, and goes through the "softest point" ^ where the ratio of the pressure to the energy 
, density takes its minimum as a function of the energy density. Therefore, it is expected that the suppression of the 
Q.^ ' collective flows is not only a signal for the existence of the QGP [^-^ but also a useful tool to determine the EOS 
near the phase transition region. 

Some experimental groups reported that the inverse slope parameters of the transverse mass spectra for the non- 
, multistrange hadrons tt, K, p, and d in central Pb + Pb collisions at 158A GeV at the CERN Super Proton Synchrotron 
■ (SPS) are parametrized by two common values, the freeze-out temperature and the transverse flow velocity [Q. This 
implies that these particles constitute the radial flow and the local thermalization at least among those particles 
is achieved in central collisions. It has been, however, an open question whether equilibration is achieved even in 
non-central collisions. In this Letter, we study the elliptic flow that is made of stable particles and resonances with a 
hydrodynamic model, and show that the final pion distribution is well described in this approach, if decay kinematics 
O , is appropriately taken into account. 

2 ' The rapidity dependence of azimuthal anisotropy for particle i is characterized by the coefficients {n — 1, 2, 3...) 
in the Fourier expansion of the azimuthal distribution of the particle S : 
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where (f) is the azimuthal angle measured from the reaction plane and Y is the rapidity. Non- vanishing vi and V2 imply 
the formation of directed and elliptic flows, respectively. Experimentally, V2{Y) is of the order of four percent around 
the midrapidity {Y = 2.92) for low transverse momentum (50 < pt < 350 MeV/c) charged pions in non-central Pb + 
Pb collisions that correspond to the impact parameter range 6.5 < b < 8.0 fm p. The value of the observed V2 has 
been claimed to be smaller than that predicted for direct pions in hydrodynamic models. 

In the following, we consider not only direct pions but also indirect pions that are from resonance decays. Thus, 
we need to strictly distinguish the matter flow before freeze-out and the observed particle flow. In this Letter we 
use a term "in-plane" when the hydrodynamic flow before freeze-out is directed preferentially to the positive and the 
negative x-axes on the transverse plane and a term "positive elliptic flow" when V2 of observed particles is positive. 
Here the x-axis is defined as the direction of impact parameter in non-central collisions. When one only considers 
particles directly emitted from freeze-out hyper-surface, "positive elliptic fiow" means that the hydrodynamic flow is 
"in-plane". However, once feeding from resonance decays is included, the above two are no longer equivalent. 

First, let us suppose that a resonance particle with mass rrifl is emitted from freeze-out hyper-surface and decays 
into two identical daughter particles with mass m in the vacuum. In the rest frame of the resonance particle the decay 
is isotropic and the daughter particles are uniformly distributed with respect to the azimuthal angle when one averages 
the spin of the resonance particles. This is, however, not the case anymore if the resonance is moving in a reference 
frame. To see this more clearly, let us assume that the resonance particle moves with velocity {VRx,VRy,VRz) = {V, 0, 0) 
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in a reference frame. If V is larger than the critical value V* = p* /E*, where p* = a/to|j/4 — ^j^d £:* = mii/2, 
the probability that the daughter particle is emitted with an angle </> from the x-axis peaks at (j)± = ±sin~^[p*(l — 
V^y/^/{mV)] and is limited in a range, 4>- < (j) < <t>+ |lO| ]. The two peaks are due to the Jacobian singularity 
in the Lorentz transformation from the resonance rest frame to the reference frame. As a result, the opening angle 
between the daughter particles tends to remain finite even if the resonance particle moves at a large velocity. This is 
the reason why the equivalence between "in-plane flow" and "positive elliptic flow" breaks down when one includes 
the decay of resonances in the final state. 

The multiplicity of pions through two body decays of resonances is given by 

Nr^^x = / -. — I Br^^x-^ I dsWR{s)—--^ ^ — . , , 2) 

J 47rp* Er J {2TTy exp[(p^u^ - n)/Tf\ =F 1 

where u'^ , p,, Tf, and ct^ are, respectively, the four-dimensional fluid velocity, the chemical potential, the freeze- 
out temperature, and the freeze-out hyper-surface, p^ is the longitudinal momentum of pion. and dR are, 
respectively, the resonance four momentum in the reference frame and the degeneracy, and — (+) is for boson (fermion) 
resonances. Br^^x is the branching ratio of the decay process and Wr is the Breit-Wigner type function, which 
takes account of the finiteness of the resonance width. For Wr^ we adopt the form used in Ref. The Jacobian 
of the Lorentz transformation from the resonance rest frame to an arbitrary reference frame J{pl,4']^ r) is defined 
by dp'^dcj}* — J{pl, 4>'j^ R)dpLd4>, where the quantities with (without) * are the ones in the resonance rest (reference) 
frame. Two typical shapes of the Jacobian as functions of the azimuthal angle in p —> tttt are shown in Fig. |l| (a). If 
Vp = {Vp, 0, 0) and Vp < V* 0.93c in this process), J{pl — 0, 4>; Vp) has a broad peak around </> = at the same 
azimuthal angle of the resonance particle as expected. However, if Vp > V*, J{pl = 0, <j); Vp) has a finite value only 
in the range < (p < 0+ and two sharp peaks appear at <j)± in addition to the original broad one as explained above. 

In order to get an idea about the effect of thermal smearing, we estimate the azimuthal distribution of pions through 
the above process, taking the following simple model: There are only two fluid elements with local fluid velocities 
[vx, Vy, Vz) = (±0.5c, 0, 0). Both elements are assumed to be thermalized at the temperature Tf = 120 MeV. The pion 
distribution from p-meson decays in each fluid element is shown in Fig. ^ (b). Here we sum up the pion distribution 
over 50 < < 350 MeV and —0.5 < Y < 0.5. The sharp peaks at ~ ±1.21 (±1.93) in the Jacobian are smeared 
by thermal motion of p-mesons, but they are still visible in the azimuthal distribution of pions from p-mesons in the 
fluid element with Vx — 0.5c (—0.5c), while the peaks at (j) — (tt) are completely washed out. The superposition of 
the two distributions leads to an azimuthal distribution with broad peaks at ±7r/2, as shown in Fig. ^ (b). In this 
simple model, pions from p-meson decay have negative W2, i.e., the elliptic flow in the final state is negative, while the 
motion of the two fluid elements is in-plane. 

We next carry out realistic hydrodynamic simulations for the space-time evolution of Pb + Pb collisions at 158^1 
GeV to see how large the effect of the Jacobian singularity is for the final state pion distribution. We first assume 
that the hot and dense nuclear matter produced in heavy- ion collisions is in local thermal equilibrium after to — 1.44 
fm/c since the two nuclei touched [p^ . Then, we describe the space-time evolution of nuclear matter after this time 
by using a relativistic hydrodynamic model without assuming the Bjorken's scaling solution |13| or the cylindrical 
symmetry along the collision axis ||l^-|l7|]. Thus, it is possible to discuss the rapidity dependence of elliptic flow 
V2(Y) for charged pions through resonance decays as well as for charged pions directly emitted from the freeze-out 
hyper-surface in this model. We use a model EOS with a first order phase transition between the QGP phase and 
the hadron phase jl6|. The QGP phase is assumed to be free gas composed of quarks with Nf — 3 and gluons. 
For the hadron phase we adopt a resonance gas model, which includes all baryons and mesons up to the mass of 2 
GeV Ip^ , together with an excluded volume correction fl^ . We use the critical temperature at zero baryon density, 
Tc{nB = 0) = 160 MeV. The two model EOS's are matched by imposing the Gibbs' condition for phase equilibrium 
on the phase boundary. The numerical results of the hydrodynamic simulation give us the momentum distribution 
of hadrons through the Cooper- Frye formula with a freeze-out energy density £'t = 60 MeV/fm^. We have used 
this formula for the direct emission of tt^, if^, and p. In addition, we have taken into account negative pions from 
the decays of p, cj, K* , and A in the final state. 

In the numerical simulation we have fixed the impact parameter at 7.2 fm. We have chosen the initial parameters in 
the hydrodynamic simulation to reproduce not only the rapidity and transverse mass distribution of negative hadrons 
in central collisions pl[ | but also the preliminary data of rapidity distribution of negative pion in non-central collisions 
p2t by the NA49 collaboration. The corresponding central energy density E^ and baryon number density ubo are, 
respectively, 3.9 GeV/fm^ and 0.46 fm"^ The energy density and the baryon number density on the transverse plane 
are assumed to be in proportion to the number of wounded nucleons with the standard Woods-Saxon distribution for 
the nuclear density. We have assumed the initial transverse flow to vanish. For details on our hydrodynamic model, 
see Ref. Q . 



2 



We first discuss the effect of tlie Jacobian singularity for the pions through p-meson decays in the reaUstic hydrody- 
namic calculation. The momentum distribution of p-mesons is free from the effects of the Jacobi function and is given 
by the second integral in Eq. i.e., the integral with respect to s and ct^. In Fig. ||, we show V2 for the p-mesons 
directly emitted from freeze-out hyper-surface and for the pions through p —>■ nn. The elliptic flow of p-mesons is 
positive and in-plane. Nevertheless, the V2 for the pions through p — > tttt almost vanishes near the midrapidity due to 
the decay kinematics. This implies that the effect of the Jacobian singularity, which we discussed above with a simple 
model, survives even in this realistic calculation. The behavior of the ?;2's of pions from K* or A is similar to this. 

Finally, we discuss the rapidity dependence of the observed pion elliptic flow obtained from our fully three- 
dimensional hydrodynamic calculation. In Fig. |^ we compare our results with the experimental data by the NA49 
Collaboration The solid line represents V2 for the total charged pions. For comparison, V2 for pions directly emit- 
ted from freeze-out hyper-surface (dashed line) and that for pions through resonance decays (dotted line) are shown 
separately. Our results were obtained by summing up pion distribution over a pr range, 50 < pt < 350 MeV/c. The 
experimental data corresponds to an impact parameter range 6.5 < 5 < 8.0 fm. The solid line is in good agreement 
with the experimental data near midrapidity ||23|] . The V2 for indirect pions from resonance decays vanishes. This 
reduces V2 for the total pions by about 26 % at midrapidity. This figure tells us that hydrodynamic description, which 
assumes the local thermal equilibrium, works well also for the expansion stage of non-central collisions at the SPS 
energy [1^. 

A few remarks on other calculations are in order here. At midrapidity, our result for the total pions is consistent 
with the previous result obtained by a two-dimensional hydrodynamic model . The authors of Ref. ||2^ assumed 
that the longitudinal expansion can be described by the Bjorken's scaling solution and numerically simulated the 
evolution of nuclear matter only in the transverse directions. Hence they could not obtain the rapidity dependence 
of elliptic flow. They also took into account resonance particles up to the mass of the A(1232), but concluded that 
resonance decays reduce the momentum anisotropy for pions by only 10-15 %. Their reduction factor corresponds to 
full pt range psf . The effect of Jacobian singularity is important in low px region. When we integrate the distribution 
over full Pt range, we obtain the result reduced by 11 % . Liu et al. p6t compared their results from a transport 
model, the Relativistic Quantum Molecular Dynamics (RQMD), with the experimental data and concluded that the 
model calculations are in reasonable agreement with experimental data. The experimental data of elliptic flow was, 
however, later updated and the agreement is not as good as before anymore. Soff et al. also obtained V2(Y) 
for pions from a microscopic transport model, the Ultrarelativistic Quantum Molecular Dynamics (UrQMD), but the 
situation is similar to the RQMD model. Recently, it was argued that the suppression of the elliptic flow is due to 
partial thermalization 29 1. However, as we have shown, this is not needed to explain the data. According to our 
result, the suppression rather implies the full thermalization of the system, including resonances. 

In summary, we have investigated the elliptic flow of pions in Pb -I- Pb non-central collisions at 158 A GeV in a 
relativistic hydrodynamic model. As sources of pions, we have considered not only direct pion emission from the 
freeze-out hyper-surface but also decays of resonance particles. Pions from resonance decays suppress the azimuthal 
anisotropy in the midrapidity region as much as 26 %. Taking this effect into account, we were able to reproduce 
the experimental data of the rapidity dependence of V2- These results lead to our conclusion that the pions and the 
resonance particles constitute thermalized in-plane elliptic flow and that the hydrodynamic picture is applicable to 
the expansion stage in the non-central collisions at the SPS energy. The Jacobian singularity effect should exist also 
in cascade calculations such as RQMD and UrQMD if the decay kinematics is appropriately taken into account. It 
will be interesting to see how important this mechanism is for the suppression in those calculations. 
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The solid and dashed hnes correspond to Vp = 0.94c and 0.9c, respectively. (b)Azimuthal distribution (in arbitrary units) 

of pions at midrapidity through decays of p-mosons from two fluid elements which are moving in opposite directions. The 
velocities of fluid elements are assumed to be {vx,Vy,Vz) = (±0.5c, 0, 0). 
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FIG. 2. Azimuthal anisotropy V2 for p-mesons directly emitted from freeze-out hyper-surface (solid line) and for pions through 
p ^ -Kw (dashed line) in non-central Pb -|- Pb collisions at 158 A GeV. 
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FIG. 3. Azimuthal anisotropy V2 for total charged pions (solid line), direct pions (dashed line), and pions from resonance 
decays (dotted line) as functions of rapidity. The experimental data was measured by NA49 in Pb -|- Pb collisions at 158 A 
GeV. See text for details. 
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